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Centaurus A (Cen A) is the nearest radio-loud AGN and is detected from radio to very high 
energy gamma-rays. Its nuclear spectral energy distribution (SED) shows a double-peak feature, 
which is well explained by the leptonic synchrotron + synchrotron self-Compton model. This model 
however cannot account for the observed high energy photons in the TeV range, which display a 
distinct component. Here we show that ~ TeV photons can be well interpreted as the ty° decay 
products from py interactions of Fermi accelerated high energy protons in the jet with the seed 
photons around the second SED peak at ~ 170 keV. Extrapolating the inferred proton spectrum to 
high energies, we find that this same model is consistent with the detection of 2 ultra-high-energy 
cosmic ray events detected by Pierre Auger Observatory from the direction of Cen A. We also 
estimate the GeV neutrino flux from the same process, and find that it is too faint to be detected 
by current high-energy neutrino detectors. 

PACS numbers: 98.70.Rz; 98.70.Sa 



/. Introduction. — Centaurus A (Cen A or NGC 5128) 
is the nearest active radio galaxy with a distance of ap- 
proximately 3.5 Mpc and redshift z = 0.002 [T]. Al- 
though by active galactic nuclei (AGN) standard its bolo- 
metric luminosity is not very large, because of its prox- 
imity to earth it is one of the best studied AGN. Opti- 
cally, Cen A is an elliptical galaxy undergoing late stages 
of a merger event with a small spiral galaxy. Sufficiently 
large amount of photometric data are available to build a 
well sampled spectral energy distribution (SED) of Cen 
A. The emission from the nucleus of Cen A has been 
observed throughout the electromagnetic spectrum from 
radio to gamma- rays [2]-[B], which shows that Cen A has 
a FR I morphology with two radio lobes, and is a non- 
blazar source with a jet inclination in the range of 15° 
to 80°. The nuclear SED shows two peaks, one in the 
far-infrared band (~ 4x 10 -2 eV) and another at around 
170 keV [7]. In the framework of the unification scheme 
of AGN, blazars and radio galaxies [8TH0] are intrinsi- 
cally the same objects, viewed at different angles with 
respect to the jet axis. The double-peak SED structure 
observed in Cen A is similar to that of blazars whose 
jets beam towards earth, suggesting that the same spec- 
tral feature is also expected from misaligned jet sources 
such as Cen A. The leading interpretation is the single- 
zone synchrotron and synchrotron self-Compton (SSC) 
model. In this scenario, the multi wavelength emission 
originates from the same region. The low-energy emis- 
sion in radio to optical wavelengths is the non-thermal 
synchrotron radiation from a population of relativistic 
electrons in the jet, while high-energy emission from X- 
rays to very high energy (VHE) gamma-rays are from the 
Compton scattering of the above seed synchrotron pho- 
tons by the same population of electrons. This model is 
found very successful in explaining the multi- wavelength 
emission from BL Lac objects and FR I galaxies such as 



NGC 1275 and M87 [i H]. Applying it to Cen A, one 
found that it is successful in explaining most of the multi 
wavelength SED data [TT] . The difficulty, on the other 
hand, is the multi- TeV emission detected by HESS dur- 
ing 2004 to 2008 [5]. Even though the HESS data alone 
can be fitted by a power law [5], and the entire 10 6 — 10 13 
eV spectrum may be still accommodated within one sin- 
gle power law model, a clear dip around 10 10 eV revealed 
by Fermi indicates an excess of TeV emission from the 
extrapolation of the Fermi data [B] (Fig[I]). So far this 
TeV spectral component is not well interpreted within 
the published leptonic models [6j [11] . 

On the other hand, Cen A has long been proposed as 
the source of very high energy cosmic rays. Recently 
Pierre Auger Collaboration reported two UHECR events 
fall within 3.1° around Cen A [12]. By assuming that 
the two events are from Cen A, the expected high energy 
neutrino event rates in detectors such as IceCube [HI U3] 
and the diffuse neutrino flux from Cen A [15j have been 
estimated. The flux of high energy cosmic rays as well 
as the accompanying expected secondary photons and 
neutrinos are calculated from hadronic models |16j . 

The astrophysical objects producing UHECRs also 
produce high energy 7-rays due to interaction of the cos- 
mic rays with the background through pp or pj interac- 
tions [Pol - f2"0"] . In this letter, we show that the multi- TeV 
7-ray emission from Cen A can be naturally interpreted 
by a hadronic model invoking pj interactions between 
Fermi-accelerated protons in the jet and the seed pho- 
tons near the SSC peak (170 keV). The same model is 
found consistent with the detection of 2 UHECR events 
from Cen A. 

II. Hadronic model of TeV gamma-rays. Observa- 
tions of variable, non-thermal high energy emission from 
AGNs imply that these sources are efficient accelerators 
of particles through shock or diffusive Fermi acceleration 
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processes. While efficient electron acceleration is lim- 
ited by high radiative losses, protons and heavy nuclei 
can reach UHE through the same acceleration mecha- 
nism. In general, these energetic charged particles (elec- 
trons and protons) have a power-law spectrum given as 
dN/dE (x E~ a , with the power index a > 2 [211 [22] . 

The dominant pj interaction is through the A- 
resonance, i.e. 

_^ ^ + _^ ( pn°, fraction 2/3 

^ ~ 1 n n + — > ne + v e v ll i> p , fraction 1/3 ' 

which has a cross section a a ~ 5 x 10 28 cm 2 . The 
charged 7r's subsequently decay to charged leptons and 
neutrinos, while neutral 7r's decay to GeV-TeV photons. 
For interactions at A-resonance, the matching condition 
is E p e' 7 ~ 0.32 (1 - cos6>) _1 GeV 2 , where E' p and e 7 are 
the proton and the background IC photon energies in the 
comoving frame of the jet, respectively. Since in the co- 
moving frame the protons collide with the IC photons 
from all directions, in our calculation we consider an av- 
erage value (1 — cos 8) ~ 1 (6 in the range of and tt). 
Denoting T as the bulk Lorentz factor of the jet, and 
T> = r _1 (l — /3cos0 o b) _1 as the Doppler factor to the 
observer (0 o b is the angle between the observer and the 
jet direction, and (3 — v/c is the dimensionless speed of 
the jet), one can re- write the matching condition as 

E p e 1 ~ 0.32 TV GeV 2 . (2) 

Here e 7 = X>e'/(1 + z) is the observed photon energy, 
while E p = TE' p /(l + z) is the energy of the proton as 
measured by an earth observer, if it could escape the 
source (instead of producing t:q photons) and reach earth 
without energy loss. This is because the proton energy 
in the rest frame of the AGN central engine (in which the 
jet is observed to move with a Lorentz factor T) is TE' p . 
An earth observer is at rest of this frame, but with an 
additional effect due to cosmic expansion. This definition 
of E p is not of significance in calculating the TeV photon 
spectrum, but is more convenient to discuss UHECRs 
(see III). In ([2]), the (1 + z) parameter has been neglected 
due to the small redshift z = 0.002 of the source. 

In the comoving frame, each pion carries ~ 0.2 of the 
proton energy. Considering each 7r° splits into two 7-rays, 
the 7r°-decay 7-ray energy in the observer frame can be 
written as E 7 = T>E' p /W = (T>/10T)E p . The matching 
condition between the 7r°-decay photon energy and the 
target photon energy is therefore 

£ 7 e 7 ~ 0.032 V 2 GeV 2 , (3) 

Modeling Cen A suggests a viewing angle Q ob between 
15° to 80° [TH I25]. The Doppler factor V is found in the 
range of 0.12 - 3.7 [Q E3 IID 121] • In the following, we 
adopt a nominal model of [6] , with the following parame- 
ters: V = 1, T = 7, the comoving blob size TZ' b = 3 x 10 15 
cm, and the comoving magnetic field strength B' = 6.2 
G, but keep parameter dependences in the formulae. This 



model fits well (the green curve in Figjl]) the bulk of ob- 
served photon flux starting from the low energy regime to 
10s of GeV. The second SED peak (SSC) is at ~ 170 
keV with the corresponding observed photon energy flux 
F 7 ( e P) = 9.0 x lO" 10 ergem- 2 s^ 1 [25 ]. It is clearly 
shown from Fig]T]that the model cannot account for the 
observations in higher energies. 

Adopting T> = 1, the target peak photon energy e v 1 = 
170 keV is matched by E* ~ 13 TeV and E* ~ 190 GeV 
for the A-resonance condition (Eq.([2| and Q). The typ- 
cial photon energy E? is in the energy range of HESS de- 
tection. The optical depth of the A-resonance process in 
the emission region can be estimated as t v1 — u'cA^-j,) 
where n 7 is the comoving photon number density in the 

jet, which is given as n'^ = r\{L 1 /T> K ) /[4-7r7\!/ 2 c(e 7 /£>)], 
with k ~ (3 — 4) (depending on whether the jet is con- 
tinuous or discrete) and 77 ~ 1. At = 170 keV, 
the observed photon luminosity is £ 7 (e 7 ) = 1-32 x 
10 42 (d/3.5 Mpc) 2 erg s" 1 . For V = 1, this gives 
n'^eP) - 1.4 x 10 6 cm" 3 (d/3.5 Mpc) 2 (^' h ,i5. 5 )" 2 ^ ,t+1 
and T pl (eP) - 2.1 x 10- 6 (d/3.5 Mpc) 2 (^' 6 ,i5. 5 ) _1 ^ _K+1 . 

The photon energy flux F 1 (EP) (effectively the 
E 2 dN/dE spectrum at E — EP) is related to the total 
proton number in the source. The total electron number 
can be constrained based on the synchrotron+SSC mod- 
eling. However, since we do not know the composition 
of the jet, especially the lepton-to-proton number ratio 
(pair multiplicity), one cannot calculate F 1 (EP) from the 
available data. We therefore derive it through fitting the 
high energy photon spectrum. 

Once Fj(EP) is fitted from the data, one can cal- 
culate the spectrum of the 7r°-decay hadronic compo- 
nent, which depends on the spectra of the protons and 
of the target photons. We assume that protons have a 
power law distribution N(E p )dE p oc E~ a dE p . Since the 
number of 7r°-decay photons at a particular energy de- 
pends on the number of protons and optical depth, i.e. 
N(E 7 ) oc N(E p )/t p1 (x N(Ep)n'Je y ) (where E 7 , E p , e 7 
satisfy the matching conditions ([2| and Q), one can cal- 
culate the 7r°-decay spectrum through the scaling 

F,(E,) <(e 7 ) (E,y a+2 

where we have used the relation E p /EP = E y /EP, 
and the power index +2 converts the photon num- 
ber spectrum to energy spectrum. We fix the pro- 
ton spectral index to a — 3.08. This index is found 
to well fit the HESS data in the high energy regime, 
and it is also the typical cosmic ray spectral index in 
the UHECR regime. At the energy EP ~ 13 TeV 
(which corresponds to the EP ~ 190 GeV peak of the 
hadronic component), the proton luminosity at ~ 13 
TeV is Lp(EP) ~ (15/2)L 7 (^)[r P7 (eP)]- 1 - 4.0 x 
10 45 erg s^id/Z.b Mpc) - 2 W b , 15.5^" l - This is smaller 
than the Eddington luminosity of the central black hole 
L E dd ~ 1.3 x 10 46 erg s _1 (M/10 8 M Q ). In order not to 
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where t v ~ 1 day is the observed variability timescale, 
which determines the size of the emission region, and the 
best-fit values of B' and V of the leptonic model parame- 
ters have been adopted. Above this energy the number of 
cosmic rays should follow an exponential decay. Bearing 
in mind the uncertainties in the viewing angle, it is pos- 
sible that the maximum proton energy can reach 57 EeV 
or even higher for a same Doppler factor T> (but with 
a larger Y). In the following, we assume that E p ^ max is 
extended to 57 EeV. 

Within our model, based on the flux at one can es- 
timate the cosmic ray flux at EP. One out of protons 
interact with the target photons to produce gamma-rays 
through 7T° decay. So the proton flux F p (E p ) at pro- 
ton energy E p is related to the high energy photon flux 
F 1 (E 1 ) at the photon energy E 1 through 



F P (E P ) = 7.5F 7 (£ 7 )[t P7 (£; p )] 



(6) 



FIG. 1. The observed spectral energy distribution E^^ 1 - 
(or F 7 ) from the Cen A core region. Colored symbol points 
are observations (with different sources marked in the figure), 
and the curves are model fits: The green curve is synchrotron 
+ SSC fit from Abdo et. al. [6], while the blue, magenta, 
and red curves are the hadronic emission component due to 
7r°-decay with Ep = 4, 13, 25 TeV, respectively. The case of 
no proton spectral break (black dashed line) is plotted for 
comparison. 



violate the energy budget constraint posted by the Ed- 
dington luminosity, it is required that the proton energy 
spectrum should break to a harder index (e.g. a ~ 2) 
at low energies. In our calculation, we introduce a break 



energy EL so that a = 2 for E p < E b and 



3.08 for 



E„ 



> E b . We numerically calculate the model spectrum. 



As can be seen from Fig[T] the hadronic model spectrum, 
along with the leptonic spectrum, can well interpret the 
observational data above GeV for a wide range of E h p 
values (from 4-25 TeV). 



For E 1 < 



E p 



the A-resonance matching condition 
(Eq.([3])) requires e 7 > e 7 = 170 keV. The drop of the 
target photon flux then results in a decreasing photon 
flux. The harder proton spectrum below E p strengthens 
the effect. 

The same applies to E y > E?, whose target photons 
have e 7 < e 7 . Since the number of protons increases 
with decreasing energy (power law distribution) , the real 
energy flux peak of the hadronic component is slightly 
smaller than EP = 190 GeV. 

III. UHECR flux. The same model can be used to 
estimate the expected UHECR flux from Cen A. The 
maximim energy to which cosmic rays can be accelerated 
is constrained by the size of the emitting region and the 
magnetic field in it. For Cen A, one has [2"9"j 



Ep^max 4 X 10 



B' \ ( U 



6.2 G J \10 5 s 



eV, 



(5) 



The factor 7.5 comes from the fact that the A-resonancc 
has 2/3 probability to decay to the pir° channel as shown 
in Eq.Q, and each pion carries 20% of the proton energy. 
At EP = 190 GeV, the best fit model flux is Ey(E») ~ 
6.6 x 10 -13 erg s" 1 cm" 2 . Since the proton energy flux 
F p oc E~ a+2 , we obtain the source proton luminosity at 
any energy E p through 



f p( e p) _ ( E. 
Fp(E P p) 



El 



-Q + 2 



(7) 



Plugging in Ep = 57 EeV, and E$ = 13 TeV, we obtain 
the UHECR flux above 57 EeV as F p (57 EeV) = 1.6 x 
10~ 13 erg s _1 cm~ 2 . 

The Pierre Auger Observatory (PAO) reported that 
there are roughly 10 UHECR events above 57 EeV con- 
centrated around the Centaurus direction, a region with 
a high AGN density [56-28 . Two of these events were 
found to fall within 3 degrees from Cen A [T^], suggest- 
ing the evidence that Cen A may be the first UHECR 
source. 

We can estimate the expected number of UHECR pro- 
ton events above 57 EeV detectable by the PAO array 
[26l|27]. Taking Cen A as a point source, the integrated 



exposure of PAO is 3 = 9000/7T km . One has to also 
consider the relative exposure w(<5) for the angle of dec- 
lination 5. For Cen A, S = 47°, and the corresponding 
value is ui(S) ~ 0.64[14|. The time duration for data col- 
lection by PAO is about 15/4 yr between 1st January 
2004 and August 2007. So the expected total number of 
UHECR proton event above 57 EeV is 

#Eve„ts=»Qa,( j )f y r = 3.7C, (8) 

where C denotes the fraction of UHECRs that can escape 
from the source region. We can see for a reasonable value 
C ~ 50%, the predicted value 1.9 matches nicely the de- 
tected 2 events from PAO. Interperting the 2 UHECRs as 
associated to Cen A implies that the intergalactic mag- 
netic fields have a strength weaker than 10 -12 G. This 



4 



is consistent with a dipole extrapolation of the galactic 
magnetic fields, but is inconsistent with a wind (r~ 2 - 
dependence) extrapolation. 

VI. High energy neutrinos. From the decay mode of 
Eq. ([!]), we can see that the fluxes of neutral and charged 
pions have the relation F(ir°) = 2F(ir + ). Since each 
neutrino shares 1/4 of the 7r + energy, while each pho- 
ton shares 1/2 of the ir° energy, the energy relationship 
between a 7-ray and a neutrino produced by protons of 
same energy satisfies E v = E 1 /2. The neutrino energy 
flux can be estimated as F v = (3/4)^+ ~ (3/8)i^o = 
(3/8)-F 7 . The maximum neutrino flux is therefore Ff, ~ 
2.5 x 10~ 13 erg s" 1 cnT 2 ~ 1.5 x KT 10 GeV s" 1 cm" 2 , 
with a typical energy Ef, = EP/2 ~ 95 GeV. This flux 
level is well below the current neutrino flux upper limit 
imposed by IceCube [3"U] . 

VII. Discussion. 

We have proposed a hadronic model to interpret both 
the TeV data and the two UHECR events detected from 



Cen A. The model requires a relatively high proton-to- 
electron luminosity ratio (of order 10 3 — 10 4 ), with the 
proton luminosity close to the Eddington luminosity of 
the black hole. On the other hand, in view of the sharp 
dependence of L p on D and the uncertainties in mod- 
eling of the lcptonic component, such a scenario can be 
validated for a choice of reasonable parameters. An alter- 
ative way of accounting for the TeV component through 
a hadronic component may be through invoking photon- 
pair cascade initiated from UHECR P7 interactions [31] 
which is beyond the scope of this Letter. 
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